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SUMMARY 


This interim report presents the investigations 
carried out by G. V. R. Rao and Associates during 
the period April 13, 1971, to June 13, 1972, under 
Contract No. NAS 2-6401 with National Aeronautics 
and Space Administration, Ames Research Center, 
Moffet Field, California. Based on reasonable 
estimates of flow conditions occurring in an axial 
fan, acoustic radiation from various noise sources 
is evaluated. Results of computations on two 
specific fans are presented, and relative signifi- 
cance of the various sources is examined. The 
information presented here constitutes a step in 
our continuing effort on understanding fan noise. 
Helpful suggestions from Mr. David Hickey and Mr. 
Brent Hodder of Large Scale Aerodynamics Branch 
are greatly appreciated. 
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pitch radius 
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velocity fluctuations 
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axial distance of the source from the fan 
center 

density 

ambient density 
phase angle 

circumferential location of the OGV wake at 
the annulus under consideration from reference 
meridial plane 
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mth component 
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directions 

fundamental component 
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root mean square value 
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I . INTRODUCTION 

The radiated noise from an axial fan is a by-product of airflow 
past the moving and stationary elements in the fan. If the 
velocity field is known, or can be appropriately estimated, 
one can evaluate the pressure fluctuations in the far field. 
Based on the work of Lighthill (1) and Curie (2) the pressure 
fluctuations caused by propagation of aerodynamically generated 
sound can be written in the form 


1 d 2 f 
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where p = acoustic pressure in the far field 
P - density 


p Q = ambient pressure 

= cartesian coordinates (i=l,2,3) with origin 
x at the source 


( 1 ) 


F i = components of the externally applied forces on 
the fluid boundary 

T^ . = components of the externally applied stresses 
J distributed within the fluid volume 

= velocity components 

and h = distance from the source. 


The subscripts i and j in the above equation indicate 
components in the three dimensional cartesian coordinate 

directions and repeated index denotes Einstein summation 
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convention. The subscript n denotes the component normal 
to the fluid boundary. The square brackets indicate that the 
quantities are evaluated at retarded times. 

The significance and interpretation of the various terms of 
Eq. (1) are discussed in ref. (3). The first term on the right 
hand side of Eq. (1) is from quadrupole sources within the 
fluid volume. Ffowcs Williams, in ref. (4) , indicated that the 
quadrupole stress density T^ ^ can be approximated by pUJUj. 

The second term in Eq. (1) represents the dipole sources located 

on the fluid boundary. These dipole sources consist of exerted 

forces F. and the momentum pU.U transferred to the fluid 
x in 

when the boundary moves with a velocity U n normal to itself. 

The last term is from monopole sources resulting from dis- 
placement of the fluid by the motion of the boundary. 

The coordinates x^^ (i=l,2,3) and distance h appearing in 
Eq. (1) can be related to R, the distance of the field point 
from fan center and \]x , the aspect angle measured from fan 
axis. It is convenient to choose the center of the rotor plane 
as fan center. As shown in enclosed Fig. (1), let us choose 
a fixed cartesian coordinate system x, ,X 2 , with its origin 
at the fan center, and plane containing the field point 
(R ,<£■). Let r, 0,£ denote the cylindrical coordinates of 
the source, which may be any one of the types described above. 

The radial distance from fan center is denoted r, whereas 9 
denotes the angle between the x 2 X 3 plane and the meridional plane 
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containing the source point. The axial distance of the source 
point from fan center is denoted £ . In the discussions 

presented in later sections 0 is considered positive when 
measured in the direction of rotor rotation. The cartesian 
coordinate system 1/2,3 with its center at the source point, 
are conveniently taken in the axial, radial and tangential 
directions as shown in the figure. By appropriate translation 
and rotation of the coordinate system we can derive the 
following relations. 


and 


x t = Rcos i J/- £ 
x z = Rsin iJ/cosQ-r 
x 3 = -Rsintfr sinQ 
h w R jl - sini/fcosQ 



( 2 ) 


In far field computations of fan noise the second order term 
in (r/R) can be ignored. 


In applying Eq. (1) to evaluate the effect the various acoustic 
sources present in an axial fan we shall adopt the following 
approach. The various elements in the fan, such as rotor blades 
stator vanes and support struts form the fluid boundaries. 

Also these elements set up a velocity field extending into the 
fluid volume. Each row of blades, or vanes, or struts can be 
replaced by a fictitious actuator disc, with the flow parameters 
and forces appropriately prescribed on them. The flow field 
set up by the blades or vanes is considered to influence the 
quadrupole strength distribution within the fluid on either 
side of the actuator disc. 



In this interim report we have considered only the quadrupole 
sources within the fluid and the externally applied forces on 
the boundaries. The influence of the remaining terms of Eq. (1) 
is presently under consideration and will be treated in a later 
report. Only such acoustic sources that contribute to the noise 
at blade passing frequency and its harmonics are examined in 
this report, since such discrete frequency noise is the prime 
origin of annoyance. In our preliminary evaluation of fan 
noise, we have ignored the influence of the duct walls. Con- 
sequently, the only boundaries of the fluid considered are 
actuator discs replacing the rotor and stator. Such omission 
is a valid approximation, if we restrict our attention only 
to the propagating modes in the duct as discussed by Tyler and 
Soffrin (5). 


-4 



II. DIPOLE SOURCES IN AN AXIAL FAN 

From the airfoil section characteristics and mean flow 
properties, the aerodynamic loads on the various blade or 
vane elements can be obtained. The forces acting at the dis- 
crete locations of these elements can be distributed over 
appropriate cross sectional planes by employing Fourier repre- 
sentation. Besides the steady state load, the rotor blades 
will ejqperience periodic fluctuating loads as they cross the 
wakes from upstream stator vanes and also random forces due 
to turbulence in the flow. The vanes in a stator located 
downstream of the rotor also experience similar steady and 
unsteady loads. In evaluating acoustic radiation from such 
dipole sources, the influence of the blades and vanes can be 
replaced by an equivalent distribution of dipole sources in 
the rotor and stator planes respectively. 


II. A. Steady state Loading on Rotor Blades 
Let us consider B number of blade elements equally spaced 
in an annulus 2-rrrdr of the fan cross section. Each blade 
element experiences a force given by 


where 


and 


A L = 


V o 


'Lr 
c _ 


I v o C Lr c r dr 


(3) 


mean flow velocity relative to the rotor blade 
rotor blade section lift coefficient 
rotor blade chord. 
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In the above representation we have ignored the drag compared 
to the lift on the blade element. Since the force A L acts 
essentially normal to the blade chord, we can write the axial 
and tangential components as 


4L, = Ah sin\ r 
and 4L 3 = A L cos\ r 

where \ r = the rotor blade stagger angle. 

The lift Al> acting on each blade element can be considered 
concentrated at the center of pressure of the blade element. 
An equal and opposite force is exerted on the fluid. As the 
blades rotate at an angular velocity Q , we obtain a fluc- 
tuating force distribution on the annulus which can be repre- 
sented by the real part of a complex Fourier series as 


where 


and 



a T i(nB0~ka o t) 

^ Be 
IT 


(4) 


n = harmonic index (=1 for blade passing frequency) 
ka 0 = n£?B 
Q = 2ttN 

N = rotor revolutions per second. 


In the above simple representation we observe that the ampli- 
tude and phase angle of each harmonic are the same, as a con- 
sequence of assuming z)L as concentrated at a point. Con- 
sideration of chord-wise distribution of A L, and its effect 
on loading harmonics is discussed in ref. (6). Amplitudes in 
higher harmonics would then fall off rapidly. 
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The term Ah in Eq. (4), as described by Eq. (3), is depen- 
dent upon the flow conditions and blade section characteristics 
occurring at the radius r under consideration. For the design 
conditions on Rao Fan 2, described in ref. (7), the radial vari- 
ation of Ah is shown in Fig. (2). The blade loading on 
LF336 fan, described in ref. (8), operating at design conditions 
is presented in Fig. (3). 

Acoustic radiation field from the forces described by Eq. (4), 
acting on the elemental annulus 2TTrdr , can be obtained by 
using the second term on the right hand side of Eq. (1). Upon 
integration over 0 , from 0 to 2-nr ^ the rms value of the 
nth harmonic of sound pressure can be written as 

J nB (krsini/0 (5) 

where J „ is the Bessel function of the first kind of order 
nB and argument krsini^ and Ah is given by Eq. (3). 

The above equation for acoustic radiation from an elemental 
annulus 2iTrdr containing B number of rotor blades, is the 
same as that derived several years ago by Gutin (9). Substi- 
tuting Eq. (3) for Ah into the right hand side of above 
Eq. (5) and integrating from r = r h to r fc , one obtains the 
radiated sound pressures from the steady state rotor blade 
loading. The results of such computations at the blade passing 
frequency on the two fans referred above are shown in Figs. (4) 
and ( 5) respectively. As can be expected from multibladed 
rotors, the theoretical estimates presented in Figs. (4) and 


Ap = 


kB 


2j2ifR 


Ah 


sinx^cos#'- cosX r 
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(5) are far below the; measured values given in refs. (7) and 
(8). Filleul in ref. (10) , also indicated the inadequacy of 
steady state rotor blade loading for predicting noise from 
rotors with large number of blades. 

II. B. Unsteady Loading on Rotor Blades 

As the rotor blades travel circumferentially, they encounter 
non-uniform flow conditions. In the following analysis we 
shall restrict our attention to the circumferentially non- 
uniform steady flow caused by upstream vanes and also the 
random fluctuations in velocity. The former gives rise to 
periodic fluctuations in blade loading, which can be treated 
as "ordered" dipoles on rotor blades in the acoustic calculations. 
On the other hand, the random velocity fluctuations present 
in the turbulent regions give rise to random dipole sources 
on the blades. Acoustic radiation from such ordered and random 
sources are examined in this section. 

II.B.l. Periodic Fluctuation in Blade Loading 
Circumferential steady non-uniformity in the flow met by 
a rotor, due to inflow distortion or the presence of upstream 
stator vanes, causes periodic fluctuations in blade loading. 

Steady state inflow distortion, met by a blade, is repeated 
in each revolution of the rotor and hence the fundamental of 
the resulting periodic blade loading is the rotor frequency. 

On the other hand, in the case of an upstream stator causing 
V number of identical wakes, the rotor blade experiences 
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periodic loading with a fundamental frequency NV. Let us 
examine the latter case in detail, since similar analysis 
also applies to the effect of inflow distortion. 


A method of estimating the lift fluctuations on the blade 
elements is discussed in ref. (11), and can be represented by 
the real part of a complex Fourier series as 


E A -i0 -iaV42t 

(4L) a e a e (6) 

a=-°o 

A 

The amplitude ( AL) and the phase angle cfr of each ath 
component depend upon frequency, wake-velocity profile, and 
the blade geometry as discussed in ref. (11). 


Due to swirl in the flow approaching the rotor, or from the 
use of non-radial upstream vanes, the wakes impinging on the 
rotor blades may have non-radial configuration. The wake- 
blade interaction in such cases does not occur simultaneously 
along the length of the blade and above Eq. (6) can be modi- 
fied to account for such effects. Based on the analysis given 
by Filotas in ref. (12), the method of evaluating the fluctu- 
ating blade loads is developed in ref. (13) and one can write 


A L 


- Yj 


-i -i(aV0t-aV# ) 


e a e 


w' 


(7) 


a= -oo 


The angle ip in the above equation, measured positive in the 
direction of rotor rotation, is the angle between the meridional 
planes containing the wake location at radius r and at hub. 
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In Rao Fan 2, upstream stators containing leaning vanes were 
employed, and the consequent values of i// at the rotor plane 
are described in ref. (7). The influence of such oblique wake- 
blade interaction on the fundamental (i.e. a = 1) component 
of the fluctuating blade loading is evaluated in ref. (14) 
for various vane-lean angles employed in Rao Fan 2. Radial 
variation of amplitude and phase angle for radial vanes and 
leaning vane configuration 3, described in ref. (7), are 
presented in Figs. (6) and (7) respectively. The large radial 
variation of phase angle shown in Fig. (7) is essentially due 
to the variation of relative velocity at the rotor in Rao Fan 2. 

Let us examine the events occurring in the elemental annulus 
2TTrdr under consideration. We have B number of blade 
elements, each experiencing load fluctuations as given by 
Eq. (7), but with a time delay of — between successive ele- 

oiN 

ments. Appropriate summation of events occurring on B number 
of blades crossing V number of wakes, as shown in ref. (14) 
would yield the following representation for the fluctuating 
force distribution on the elemental annulus. 

oo oo 

£ £ <-> 

n= 1 a= -oo 

where ka 0 = 2'rmBN 

n = harmonic index ( =1 for blade passing frequency) 

m = nB-aV, but takes on values between -°o and +°o (9) 


i(m0- 6 ) 

r ma 


-i(ka 0 t-aVi^) 


( 8 ) 
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Above Eq. (8) is similar to Eq. (4) derived earlier for the 
steady state loading on rotor blades. However, we observe 
that the nth harmonic of the above impressed force contains 
contributions from all the ath components of the fluctuating 
load on the rotor blades. Taking the axial and tangential 
components of the externally applied force described by Eq. (8) 
we obtain 



dF, 

= dF 

sin\ r 

and 

dF 3 

= dF 

cosA. r 

where 

^r 

= the 

rotor blade stagger angle 


( 10 ) 


The acoustic radiation from the dipole sources equivalent to 
the fluctuating forces described above, can be examined using 
the second term of Eq. (1). It is shown in ref. (15) that the 
acoustic pressure at the field point (R,tf0 is then given by 


^ p n 4 ttR 



a=-°o 


/ mcos\ 

( sin\ cos ilf 

V r 


kr 


(-i) m+1 


i(aV ip - 0 ) 

X (it) & e W ma J^krsintf' - ) 


/\ 


( 11 ) 


with the values of m given by 
m = nB - aV. 


Since the exponents of 0^ and 0 depend upon radius r 
and the summation extends over several values of m and a, 
we shall abstain from expressing the rms value of pressure 
described above. The above equation agrees with the represen- 
tation given by Lowson (16). Integrating right hand side of 


-11 



Eq. (11) over r^ < r < r fc , we obtain the contribution to fan 
noise from the dipole sources due to wake-blade interaction 
forces. The results of such computations on Rao Fan 2 are 
given in ref. (14). The directivity of sound pressure level 
for the typical vane configurations is presented in Fig. (8). 
The contribution to fan noise from periodic blade loading is 
several orders of magnitude larger than that from steady 
state blade loading, shown in Fig. (4). This is a consequence 
of the Bessel functions of order m in Eq. (11) and NB 
in Eq. (5). The significance of wake-blade interaction in 
fan noise generation depends strongly on blade-vane number 
ratio. Even though the use of leaning vanes decreases side- 
line noise as indicated in Fig. (8) and confirmed in measure- 
ments reported in ref. (7), the periodic blade loading con- 
sidered here is inadequate to predict fan noise at on-axis 
points. 

The effect of inflow distortion on blade loading and resulting 
noise can be obtained by substituting V = 1 in Eq. (6), and 
the amplitude (dL) a depends on Fourier series representing 
the distortion, since inflow distortion is more diffused 
than the narrow wakes of stator vanes, as the index a 
increases the amplitude (4L) decreases rapidly, a double 
Fourier series in r and 0, to represent the flow distortion 
yields the corresponding value of (dL) as a function of radius 

Before closing the discussion on periodic fluctuations in blade 
loading, it is interesting to note the influence of the wakes 
from support struts in LF336 fan impinging upon the downstream 
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rotor blades. There is a major strut spanning the fan annulus 
and also aMninor strut perpendicular to it. The noise generated 
by the presence of the minor strut can be ignored relative to 
that due to the major strut, as discussed in ref. (17). The 
major strut spanning the inlet can be considered as two radial 
vanes upstream of the rotor and one can employ Eq. (9) to evaluate 
the radiated acoustic field. In using Eq. (9) for acoustic 
computations on Rao Fan 2, it is found that a single value of 
a = 1 is sufficient. Higher values of a, and associated 
higher orders m of the Bessel function led to negligible changes 
in the computed results. However, in the case of LF336 fan with 
B = 42 and V = 2 ( corresponding to one strut across the fan 
annulus) , all values of a ranging from 1 to 21 had to be 
included and such computations were reported in ref. (17). 

To illsutrate the radial variation of the fluctuating loads on 

the rotor blades in LF336 caused by the major strut the ampli- 
A 

tudes ( 4L) and phase angles 0 are presented in Figs. (9) 

and (10) respectively for values of a = 1 and 21. To 

illustrate the dependence of each harmonic component of A L 

on the index a, the vector 

A -i 0 a 
UL) a = UL) a e a 

at the hub location of the blade for all the values of index 
a considered is shown in Fig. (11). The variation of the 
phase angle <f> with the index a is primarily due to the 

a 

frequency of the ath component of the fluctuations in the flow 
conditions met by the rotor blade as discussed in ref. (11). 
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Acoustic computations based on the influence of the major 
strut in LF336, are presented in Fig. (12). The estimated 
sound pressure level at the field point on fan axis is within 
reasonable agreement with measurements reported in ref. (8). 
However, these estimates of the effect of major strut fell far 
below the levels measured on the side line, indicating that 
radiation from other acoustic sources needs to be considered. 


II. B. 2. Random Fluctuations in Rotor Blade Load 
The flow conditions experienced by a rotor blade element, due 
to turbulence in the inflow to the fan, may contain randomly 
fluctuating component v Q . In the following we consider the 
random loading on the blades and acoustic radiation therefrom. 
Based on the discussion presented in ref. (18) the power spec- 
trum of the turbulent fluctuations of v , component normal to 

u n 

blade chord, can be represented as 


4V 0 <“> = v o 2 


£ l+3(^^v 0 ) 2 


TTV 0 


|i+(^o//v 0 ) 2 j 2 


( 12 ) 


where 


= mean flow relative to blade 


and £ - length scale of turbulence. 

The fluctuating cross-component of velocity would give rise 
to fluctuating lift force on the blade. Based on an approxi- 
mation suggested by Liepmann (19) for the lift response function, 
we can write the power spectrum of the fluctuating lift force 


per unit span of a blade as 

* 

<£> ItO ) = ( 7T P V 2 C f 

' ' 1T y o w o c r l+Tr c r ^'/v 0 


(13) 
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rotor blade chord 


where c = 
r 

and 3> = the spectrum described in Eq. (12). 

v o„ 


If the random fluctuation occurring on one blade element has 
no correlation with that on other elements on the blade and 
also with that on other blade elements in the fan annulus, the 
acoustic radiation from such random forces will give rise to 
broad band noise. However, if the length scale of the turbulence 
is such that an eddy is crossed by several blades before it 
loses its identity, we can observe sharp peaks in the noise 
spectrum at the blade passing frequency and its harmonics. 

This number of possible crossings is given by 


0 = largest ! integer smaller than BN- tt — ( 14 ) 

°' 1 

where = length scale (axial direction) for v 0 fluctuations. 

and U 3 = axial velocity considered here to be the same as 
1 the velocity with which the eddy is convected 
downstream. 

Depending upon the length scale in the tangential direction, 
it is possible for more than one blade to experience the same 
conditions simultaneously. The number of such blades is given 

by 


0 , = lowest integer larger than 


A 


(15) 


where = length scale (tangential direction) for v ft fluctuations 
and d = blade spacing 

n 

To facilitate our computation, we made a simplifying assumption 
that the velocities occurring over a span length -£ 2 (the 
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turbulence scale in the radial direction) are the same as those 
occurring at the center of the eddy. 

Taking into account axial components of the random lift forces 
on all blade elements in the fan annulus from r = r^ to r fc , 
it is shown in ref. (20) that the power spectrum of radiated 
acoustic pressure at field point (R,U0 can be written as 

2 

CO 




<D p ((j). = 


B 

16tt 2 R 2 a 2 . 


co 

A 2: 


m= 


— —CO 


J 2 (~ rsin iU) f 2 ( X)$> (m Q- w) ( dr 
m a n f 


(16) 


where f ( ) = 


l-cos2TTmg 3 /B ] ma 0 cos\ 

1-COC2WB / ’ fU) = 

and $ (mi3-6J) = function described by Eq. (13) at 
u (jj‘~ (mfl-u), 



1 -cos 2 tt g t cj /B S2 
1-cos 2 it oj /B j2 




Evaluation of the right hand side of above Eq. (16) becomes 
simple for axis points, since we set ^ = 0. The results of 
power spectrum of acoustic pressures thus computed on Rao Fan 2 
and LF336 are presented in Figs. (13) and (14). These preliminary 
computations show the effect of increased length scale of 
turbulence as reported by Mani ( 21) . The band width at blade 
passing frequency resulting from our computations on the axis 
appears to be narrower than that reported in ref. (21) for the 
total acoustic power from the fan. To illustrate the directivity 
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of discrete tone noise at the blade passing frequency, we 
evaluated the right hand side of above Eq. (16) at various 
aspect angles and the results for Rao Fan 2 are shown in 
Fig. (15). 

II. C. Unsteady Loading on Outlet Guide Vanes 
Since the loactions of the stator vanes are stationary, only 
the unsteady loads on the vanes need be considered. These 
unsteady loads can originate from periodic interaction with 
wakes from upstream blades or from random velocity fluctuations 
in the turbulent regions. 


II.G.l. Periodic Fluctuations in Vane Loading 
Outlet guide vanes positioned downstream of a rotor are sub- 
jected to unsteady flow conditions due to periodic impingement 
of rotor wakes. As in our earlier discussions, let us direct 
our attention to an annulus 2trrdr containing V number of 
stator vanes downstream of B number of blades. The primary 
influence of the velocity deficit in each wake is to cause 

\ 

fluctuations in the incidence angle of mean flow on the vane 
element. Since a vane element is subjected to B number of 
wakes during each rotor revolution, these fluctuations will 
have a fundamental frequency BN, and the resulting fluctuating 
lift on the vane element can be expressed as the real part of the 
complex quantity 



: /\ -i0 -inBi2t 

UL) n e n e 


( 17 ) 
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where the index n = 1 corresponds to the blade passing fre- 
quency. Above Eq. (17) is similar to Eq. (6) presented earlier 
for blade loading harmonics, except that the fundamental fre- 
quency here is BN in place of VN employed in Eq. (6). 


Due to the swirl in the flow downstream of the rotor, and due 
to non- radial disposition of the blade trailing edge, the wakes 
from the rotor blade would occupy non-radial positions as they 
reach the stator plane. Consequently, wake impingement on a 
vane does not occur simultaneously along its length. To increase 
the obliquity of wake-vane interaction, the vanes in the down- 
streamiistator can be made to lean in the direction of rotor 
rotation. Modification of above Eq. (17) to take into account 
such oblique wake-vane interaction is discussed in ref. (15), 
and the fluctuating lift on the vane element can be described 
as 


AL = 


oo A 

x A 

E UL) n 



-inB£?t 

e 


-inB( ip -ip ) 
' w ’v 
e 


(18) 


The angles 0 w and tp^ in the above equation denote respec- 
tively the angular locations of the wake and the vane at radius 
r measured from the meridional plane containing the vane 
location at the hub. In the notation employed here both ip < 

and ip^. are measured positive in the direction of rotor rotation. 

A 

The method of evaluating the amplitude ( /4L) and phase angle <p 

n n 

of each nth harmonic component in above Eq. (18) is discussed 
in detail in ref. (13), and the obliquity of wake impingement 
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plays an important role in determining the fluctuating lift 
load on the vanes. When the outlet guide vanes are radial, we 
can set 1/^=0 in Eq. (18) and obtain a representation 
similar to Eq. (7) derived earlier. The exponent containing 
\b has opposite sign in Eq. (18) compared to earlier Eq. (7), 
since the rotating wakes impinge on stationary vanes in the 
present case. 

Computations were carried out on the 45 radial vane configur- 
ation of LF336 for the amplitude and phase angle of the funda- 
mental component (n=l) of the fluctuating load, represented by 
Eq. (18) above, and the results are presented in Figs. (5 ) 
and (16). As the rotor-stator gap is increased to twice the 
rotor chord length, the severity of wake-velocity deficit is 
reduced and the obliquity of wake impingement increases. The 
combined effect is that the magnitude of the fluctuating load 
on the vanes is considerably reduced. With the radial vanes 
we have = 0, but ^ increases with rotor-stator gap. 

Calculations on various OGV configurations on LF336 are pre- 
sented in ref. (14). 

Let us assume that all the vane elements within the elemental 
annulus 2iTrdr are identical and the flow conditions in each 
of the rotor blade wakes are the same. Even though each vane 
element experiences lift fluctuations described by Eq. (18), 

the events on successive vanes occur with a time delay ^ , 
which is the time taken by a wake to travel from one vane to the 
next. We can consider the force AL described by Eq. (18) 




- 19 - 



as concentrated at each vane location, but the fluctuations 
containing a delay term consistent with its location. 

Acoustic radiation from an array of such concentrated forces 
(dipole sources) as described above is examined in ref. (15). 
The nth harmonic component of sound pressure at the field 
point (R,tfO due to such sources in an elemental annulus 
2iTrdr is shown to be 

- ika » (t - R / a o> -i0 n iKa 0 r 

^Pn ■ 45r UL) n e 2s e 

>1 

-ikrcos© . sin^ I 

-cos\ s sinQj sintfO e J j (19) 

where 9 j = 2iT(j-l)/V, and r = 0j/2iTN 

In the above equation, we observe that the various functions 
appearing on the right hand side depend upon the radial loca- 
tion of the sources considered. After summing over the V 
number of vanes and integrating from r = r^ to r fc we can 
evaluate the rms value of acoustic pressure at the field 
point (R,^). The results of such computations on LF336 
operating at design conditions is presented in Fig. (17) as 
SPLdB at 150 ft. distance from fan center. By integrating the 
acoustic intensity over i^= 0 to 180, we obtain the total 

power of radiated sound. The power levels thus estimated, at 

the blade passing frequency are presented in ref. (14) and 
show good agreement with those calculated from measured data. 
However, the directivity of radiated acoustic field as presented 
in Fig. (17) is far different from the measurements reported in 


|( sin\ s costf/- 
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ref. (8). Also, the ordered dipole sources arising from the 
wake-vane interaction considered here yield zero sound pressure 
along the fan axis. Consequently, it appears that the noise 
sources considered here, by themselves are insufficient to 
predict fan noise. 


Before closing the discussion on periodic fluctuations in vane 
loading, it is interesting to note that the concentrated loads 
at the vane locations can be replaced by a continuous distri- 
bution with appropriate Fourier coefficients. Such represen- 
tation was derived in ref. (22), and the wake-vane interaction 
effect is equivalent to the following force distribution exerted 
on the elemental annulus 2iTrdr at the stator loaction. 

n= 1 m=— 



X 


-i(m-nB) </> v J 


( 20 ) 


where ka 0 = 2'TTNBn 

m = nB - aV 

and a = all integers from -°o to +oo including zero. 
In terms of the index a, one can rewrite the above equation 
as 


oo oo 

= EE 

n=l a=-oo 
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and observe the similarity with previous Eq. ( 8) derived for 

force distribution in rotor plane. The axial and tangential 

components of the force described above are respectively 

dF, = dF sin\ 

1 s 

and dF 3 = -dF cos\ s (22) 

where \ = outlet guide vane stagger angle 

b 

The tangential component here has a sign opposite to that of 
the corresponding term in Eq. (10), since the torque on the 
OGV stator is opposite to that on the rotor. Using these 
externally impressed forces in the second term of Eq. (1), the 
radiated acoustic field can be evaluated. As shown in ref. (22), 
the nth harmonic of the acoustic pressure at the field point 
contributed by the forces on elemental annulus 2iTrdr of the 
stator plane can be written as 


_ _kV_ „ 
^ ^n 4 ttR 


-ika Q (t-r/a 0 ) 


^ | (BinA. s cos\p-+ ££ cos\ s )(-i) m+1 

a— — '*■ 


a=-°o 
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X UL) n e 
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(23) 


where m = nB - aV. 


One observes the similarity of this equation with Eq. (11) 
derived earlier for fluctuating loads on rotor blades. The 
term containing cos\ has a change in sign, because of the 
opposite sense of torque on the OGV stator relative to that 
on the rotor. The above equation agrees with the representation 
given by Lowson in ref. (16), except that the non-radial orien- 
tations of the wakes and the vanes are included here. Since 
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all the terms on the right hand side of Eq. (23) are dependent 
upon r , one completes the integration over r from r = to 
r fc before evaluating the rms value of the acoustic pressure 
at the field point. Typical computations were carried out on 
LF336 using above Eq. (23), The difference between the results 
obtained by using Eq. (23) and the earlier Eq. ( 19 ) is not 
discernible. 

II. C. 2. Random Fluctuations in Stator Vane Load 
Each stator vane element experiencing turbulent velocity fluc- 
tuations , whether in the inflow or the wakes of upstream elements, 
would suffer random lift forces. However, since the vanes 
occupy fixed locations, the acoustic radiation from such random 
dipoles would have the same spectrum as the turbulence energy. 

Our present effort is in evaluating fan noise at the blade 
passing frequency and its harmonics. Hence, the investigation 
of the influence of turbulence on stator vanes is beyond the 
scope of this report. 
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III. QUADRUPOLE SOURCES IN AN AXIAL FAN 


The influence of quadrupole sources can be examined from the 
first term on the right hand side of Eq. (1). By carrying 
out the differentiation with respect to , Xj taking into 

account the retarded times, Lighthill (23) has shown that 
the pressure in the radiated field can be written as 


p-p = 

° 4Tta 


k/ 5 


dV 


(24) 


The integration on the right hand side is to be carried out over 


the volume where T„ non-zero, and the square bracket 
indicates that the quantity is evaluated at retarded time. 


As noted earlier, we can approximate 

T ij = ^i D j < 25 > 

where U is the velocity in the region with components 
i,j = 1,2,3. Before we can utilize the above Eqs. (24) and 
(25), we should examine the flow field in the fan annulus, i.e. 
obtain an appropriate description of velocity U. As mentioned 
earlier, the velocity field can be postulated in terms of the 
flow surrounding the blades, vanes or any other bodies immersed 
in the flow. 


III. A. Description of Flow Field 

Let us represent the velocity field as 

U = U + U + U + U + U + v 
o pr ps wr ws 


(26) 


where U Q denotes the steady state mean flow which may not 
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to be uniform in the fan annulus. The subscript p denotes 
the velocity field set up by the potential flow around the 
immersed bodies and the subscript w denotes that caused by 
the velocity deficit in the viscous wakes. The second subscript 
r or s refers to the contributions from rotor blades or 
stator vanes respectively. The last term v denotes random 
fluctuations, of velocity. 

For the sake of simplicity we shall ignore any distortion 
the mean flow at the fan inlet. Also the random fluctuations 
v are thought of in terms of only the turbulent fluctuations 
and not those arising from blade-stall or other sources. 

Each component, except v , in the above equation is well 
correlated i.e. the value at any one point is definable in terms 
of its value at any other point. On the other hand, the tur- 
bulent fluctuations v may have correlation only in a small 
region. 

III.A.l. Steady State Mean Flow 

Except in the region of the inlet bell— mouth, the radial com- 
ponent Uq 2 can be ignored. The tangential component U 0;j 
is due to swirl behind the inlet guide vanes or behind the 
rotor in an OGV fan. Its magnitude will be much less than the 

axial component U 0 ^ . 
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III. A. 2. Perturbations from Potential Flow 

Let us consider an elemental annulus 2-rTrdr of the fan contain- 
ing B number of blades or V number of vanes. The potential 
flow around these blade or vane elements can be obtained in 
terms of two-dimensional flow past a cascade of infinite number 
of airfoils. The cascade geometry and the uniform velocity 
through the cascade would be those corresponding to the condi- 
tions occurring at the annulus under consideration. Each air- 
foil in the cascade can be replaced by suitable distribution of 
vortices, sources and sinks to evaluate the flow field. Such 
analysis was presented in refs. (24) and (25). It is found that 
a description of the flow field upstream of the leading edges 
is simple. In between the leading and trailing edges abrupt 
changes in magnitude and phase angle occur. The perturbations 
from the mean flow are relatively small in the region downstream 
of the trailing edtes. 


Based on refs. (24) and (25) we can write the following expression 
for the axial component of the nth harmonic of velocity in 
the region upstream of rotor leading edge. 


A 


U pr = 
* in 

U 

pr 
* in 

e 

where 

A 

U 

pr ,n 

and 


and 


(27) 


0 r = the amplitude and phase angle of 
p in the component at blade leading edge 

x* = distance measured upstream from rotor 
leading edge 

d^ = blade spacing. 
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depend upon 


The amplitude U and the phase angle 0 

in pr m 

the blade profile, blade loading, solidity and stagger angle. 

The above simple representation for U is applicable only 

pr in 

upstream of leading edge since 0 does not depend upon x' 

pr m 

in this region. Computations at pitch radius on Rao Fan 2 ' 
indicate that for the corresponding tangential component 


b = ?' + 
pr m 2 


Consequently, we can write 


iTr/2 


(28) 


Within the fan annulus the radial components will be relatively 
small and can be ignored in evaluating the quadrupole radiation. 

In the case of a stator, the flow around the stationary vanes 
gives rise to velocity perturbations which are periodic in the 
circumferential direction. The fundamental period would be the 
vane spacing and the nth component of the perturbations related 
to the stator can be written as 


a PS fin -2mTx'/d inve 

U = U e : n e s e 

ps ps B 

r in ^ in 


U ~ U e 
pS 3 n ps m 


irr/2 


(29) 


(30) 


where x' = distance measured from stator leading edge. 
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We observe the similarity between Eq. (27) derived earlier for 
the rotor related potential flow, and the above Eq. (29), 
since the only difference lies in translation of coordinate 
system with an angular velocity Q . As in the case of rotor 
related flow field, here also we can ignore the components in 
the radial direction. 

The method of evaluating flow perturbations arising from 
potential flow around the rotor blades and stator vanes, is 
discussed in detail in refs. (24) and (25). Non-radial orien- 
tation for rotor blades is seldom employed and Eq. (27) is 
suitable for all radial locations. We can substitute 0 + 0 v 
in place of 0 in Eq. (29) to include the effect of non-radial 
vanes, with t p described as a function of radius r. 

III. A. 3. Perturbations Due to Velocity Deficit in the Wakes 
In the region downstream of stator vanes or support struts, 
there are present wakes with a velocity deficit relative to 
the mean flow through the fan annulus. The representation of 
flow downstream of a row of support struts is similar to that 
in the case of stator vanes, except that there are fewer number 
of wakes and the struts produce no swirl in the mean flow 
through the fan. 

Let us consider an annular region 2iTrdr downstream of a row 
of V number of stator vanes. This region will contain V 
number of wakes equally spaced along the circumference. The 
circumferential location of each wake at the annulus under 
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consideration will not be the same as that of the corresponding 
vane element in the stator , if there is any swirl present in 
the mean flow through the annulus. Let s , measured from a 

T W 

reference meridional plane, denote such circumf erential location 
of a wake at the annulus under consideration. Let ( 5V)° denote 
the velocity defect at the wake center line and Y denote the 
wake half width. Both ( 6v)° and Y can be estimated from the 
vane section characteristics, mean flow velocity at the 
stator, and the downstream distance from the stator. One can 
assume that across each wake there exists a velocity profile 
similar to that occurring in the wake of an isolated airfoil 
section. It is shown in ref. (11) that the integrated effect 
of V number of such wake can be represented by the real part 
of a complex Fourier series 


(5V) S 

where 



v 



Jff iaV <9-<0 -< a «A s ) 2 

„ e e 


Vrrd 

= -y— cos z/ s 

= circumferential spacing of the wakes, same as 
that of the stator vanes 

= angle between the wake center line and normal 
to the elemental annulus 


(31) 


We observe that (6v) on the right hand side of the above 
equation is measured in the same direction as the velocity v , 
Hence, the components of the velocity perturbations caused by 
stator wakes can be written as 



( 32 ) 


U 


ws. 


= (6V) cos v 


U 


ws. 


= ( 6V) sin v 


From equations (31) and (32) we observe that these velocity- 

components depend upon the axial , radial and circumferential 

locations of the point, where they are to be evaluated. Since 

we have ignored radial components in our analysis, we can set 

U s 0. 

WS 2 

The analysis of flow field downstream of B number of blades 
will be similar to the above, if we consider the coordinate 
system attached to the rotor. However, in. fixed coordinate 
system we obtain 


<6V) r = 2<6V>; ^ 


r 2 

inB(e- ^ w ) -(mTA r ) -i*a 0 t 
e e e 


(33) 


n=l 


where ka Q = 2irnBN 


V'fTdj. 

K r = — — cos 

d r = circumferential spacing of the wakes, same as 
that of the rotor blades 


v. 


= angle between the wake center line (in relative 


flow) and normal to the elemental annulus 
Similar to Eq. (32), the velocity perturbations caused by rotor 
wake can be written as 


U 


wr, 


U 


wr. 


U 


wr, 


= ( 6V) cos v. 


= 0 


= ( 6V) sin v 


(34) 


with (6V) r as defined by above Eq. (33). We observe that the 
index n = 1 corresponds to blade passing frequency. 
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In our analysis here, both 0 and are considered positive 

when measured in the direction of rotor rotation. The equations 
derived in this section are for velocity perturbations at any 
given point (r,0,x) in the fan annulus and as such the respective 
parameters at the location should be used. 

III. A. 4. Random Velocity Fluctuations in Turbulent Regions 
The velocity perturbations discussed in previous subsections 2 
and 3 are well-ordered in the sense that the velocity at one 
point of the flow region has a definite space-time relation 
with the velocity at another point of the flow region. Any 
complete description of flow field cannot preclude the presence 
of random turbulent fluctuations. These random fluctuations 
can be present in the flow entering the fan, and also can be 
generated by various elements of the fan such as struts, rotor 
blades, stator vanes. Because of their very nature, we can 
only refer to these random fluctuations in terms of statistical 
quantities such as mean square values, probability density 
functions, correlation functions and power spectral density 
functions. In the study of fan noise problems, another impor- 
tant quantity is the correlation length, £ , which is defined 
as the area under the non-dimensional correlation function and 
thus represents some average size of the eddy within which two 
points in the flow field will have velocity fluctuations that 
can be considered as coherent or related. 
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The spectral distribution and the correlation length are not 
independent of each other. This is obvious for the reason 
that if low frequency fluctuations are dominant, the extent 
to which two points in the flow field are related will be 
larger. The way they are related depends on the particular 
model chosen to describe the turbulent field, as discussed 
in Appendix II of ref. (18). 


For inlet turbulence, the assumption of homogeneity and isotropy 
is usually made. In other words, the mean square values of 
the fluctuating velocities do not depend on special location 
nor the coordinate system chosen. If a "frozen pattern" 
convected with velocity U c is chosen, as discussed in Ref. (19), 
the spectra of the velocity components can be written in terms 
of longitudinal correlation length scale and mean square 
value v* of the inlet turbulence 



We note that v? = v? for isotropic turbulence. The above 

o, o 2 


spectral density functions are derived from an assumed expon- 
ential form for the longitudinal space correlation, and these 
were used in the dipole calculation. However for the quad- 
rupole calculations, the space-time correlation function of 
the turbulent field is required. Based on the same model used 


before, it has been found that the space-time correlation for 

-32- 



the longitudinal velocity component is too complicated for the 
quadrupole evaluation. Thus a slightly different model has 
to be adopted. Without assuming isotropy, the space-time correl 
ation function of the "frozen pattern" conr/ecfced with velocity 
U in the 1-direction can be written as 


c 

*r u c T 





R(* f r) = V , 2 e 

A 


-4 


A \ 


Vhere £ is 'the separation vector with components £ £ and £ 3 . 


If one is interested in that part of the noise generated by 
the true temporal fluctuation of the turbulent field, an addi 
tional time function has to be incorporated as explained in 
ref. (18). 


For wake turbulence behind struts, blades or vanes, some 
phenomenological approach can be used to estimate the intensity 
level in terms of the mean flow and some properties of the air- 
foil. Ref. (26) shows that the root-mean-square value of tur- 
bulence can be written as 

v = 0.2 (6V)° (38) 

and the length scale can be assumed as the mixing length leading 
to 

1 = 0.2Y (39) 

The velocity deficit (6 V)° at wake center and wake half-width 
Y are discussed in the preceeding subsection. 



III.B. Ordered Quadrupoles Due to Periodic Velocity Fluctuations 
Let us examine the consequence of the periodic components of 
the velocity fluctuations described in the preceeding section. 


These velocity fluctuations are well ordered, i.e. the fluctu- 
ations at any two different points in the flow region have a 
definite phase relationship between them. Let us consider a 
point (r,0,£) in the fan annulus and an elemental volume 
rdr dQ d£ within which the velocities can be assumed to be the , 
same as at the point considered. Restricting our attention 
only to the axial components of the fundamental of the periodic 
fluctuations of velocity, we substitute 


U = U + U 

o pr 


m 


+ u + u + u 
P s ,(.) wr .(.) ws '(0 


(40) 


into the right hand side of Eq. (24). The second subscript 
(1) in the above equation indicates the fundamental of the 
spatial or temporal perturbations. From the fundamental com- 
ponents of the velocity fluctuations, we can evaluate acoustic 
radiation at the blade passing frequency. The analysis can be 
easily extended to higher harmonics. 


III.B.l. The Effect of Rotor Related Flow Field 
In the region upstream of the rotor we have the velocity fluc- 
tuations U and in the downstream region we have the velocity 

P r i 


fluctuations U 


wr, 


From Eq. (27) we can write the following 


expression 
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( 41 ) 


U = 0 cos(B0-ka o t ) 

P r i(i) P r i(i) . . 

A A A ) -2TTx /d r 

with U = /U > e 

p '"> « P,W ke. 

where x‘ = upstream distance from blade leading edge, and 
the subscript _^.e. denotes the amplitude at 
the blade leading edge. 


Up to a downstream axial distance of two chord lengths velocity 
fluctuations caused by the wake and described in Eqs. (33) and 
(34) can be approximated as 
U 


= U cos(BQ-kat) 

wr,^ wr,M 0 


>(0 


(42) 


where the amplitude 
= 4,84 cf V 
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'Dr 


wr 
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r o 


0 . 6 + 
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0.15 + 
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c cosy 
r r 


c r cosy r 


with x" = downstream axial distance from blade trailing edge 
and V Q = velocity relative to rotor blades. 


Acoustic pressure due to radiation from the "ordered" quadrupoles 
in the upstream region of the blades can be written as 


r, 2 t r 

t oo 


P H) - p » = i^5 Iff cosV 


r h ° ° 


* f 
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2pU 


V,J 


rdrdQdx 1 


(43) 


where U is described by Eq. (41) and the square bracket 

pr i(0 

indicates that the function is evaluated at retarded time. Com- 
pleting the integration over 9 we obtain the rms value of 
the contribution from the elemental annulus 2'n‘rdr extending 
from x‘ = 0 to oo as 
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dx* 


(44) 


/ ’ U °, cosV J B (krsin,;,)rar / V lfl . 

o ^ 

where U is defined as a function of x* by Eq. (41). 
P r i(i) 

Similarly the acoustic pressure due to radiation from the 
•ordered* quadrupoles in the downstream region of the rotor 
can be written as 



1 

4rr a 2 0 R 


2tt 2c r 





rdrdedx" 


(45) 


where U is described by Eq. (42). 

wr '« 

The rms value of the contribution from the elemental annulus 
2iTrdr extending from x" = 0 to 2c r can be written as 

2c r 

4P(„ =575^0, J B ikrsin ' 1 ' )lAl / V |(l) dx “ (46> 

o 

let us consider only the axial component of the steady state 
blade loading and rewrite Eq. ( 5 ) as 


4P W ■= 2 SS 72 Ah si,u r V*” 11 "" > 


(47) 


We observe certain amount of similarity in the expressions 
given in Eqs. (44), (46) and (47) for the contributions from the 
various sources. For example, comparing Eq. (44) with Eq. (47) 
we can write 
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00 


(48) 


4P(,) from Eq. ( 44 ) k/>U 0] cos ijs 


4P ()) from Eq.(4 7) 


B Ah sin\ 


0 


Pr, 


dx' 


•r X r "iO) 

Since Ah sin^ r is the axial force per blade element in the 
annulus under consideration , let us define a thrust coefficient 
C T by- the equation 


Ah sin\ r = | Vq C t c r dr (49) 

The total strength of quadrupoles present in the volume con- 
sidered divided by the number of blades gives us quadrupole 
strength per blade and let us define a coefficient by 

the following equation. 

00 

i "r dr = f rdr fJ P r v , < 5 °> 

Subscript (44) on denotes that the evaluation is in terms 

of quadrupoles employed in deriving Eq. (44). 


Substituting Eqs 

ip ()) from Eq.(4 4) 
AP ()) from Eq .( 4 7 ) 


(49) and (50) into Eq. (48) we obtain 
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O, 


C V 
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(51) 


The comparison made above was also observed by Ffowcs Williams 
and Hawkings (4). The importance of the ordered quadrupoles 
increases with the number of blades and the tangential velocity 
°f tbe blades. By comparing Eq. (46) with Eq. (44) we can write 


4P (1) from Eg^46) 
AP (l) from Eq.( 44 ) 




(52) 
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Considering the magnitude of blade section drag coefficient, 
we can conclude that "ordered" quadrupoles present in the 
region downstream of a rotor are of the same importance as 
those present in the upstream region. 


III.B.2. Effect of Stator Related Flow Field 

The presence of stator vanes give rise to perturbations to the 
mean flow through the fan. As discussed in an earlier subsec- 
tion, these perturbations have space-wise dependence, but are 
of steady state nature. Consequently, the stator related flow 
perturbations in themselves produce no noise. However, the 
presence of periodic fluctuations due to rotor flow field 
coupled with the space-wise perturbations due to stator can 
generate noise at the blade passing frequency and its harmonics. 


Let us examine the case of outlet guide vanes for example. 
Using the analyses developed in the preceding section, in the 
region upstream of the stator the fundamental component of 
axial velocity perturbations can be written as 
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P s >(0 


= u 


ps 


COSV 0 
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where V = number of outlet guide vanes, and 
a f a 1 -2TTx'/d, 

u = Ju 


(53) 
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e . 

and x' = upstream axial distance from stator leading edge. 


In the region downstream of the stator up to two chords length, 
the axial component of the fundamental of the perturbations can 
be approximated by 
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0.15 + 


X" 


C COSJ/ 
s s 


V o 


= velocity relative to stator vanes 


and 


x" = downstream distance from stator trailing edge, 
s 


The fundamental of the axial velocity fluctuations caused by 
the rotor wake are already described earlier in Eq. (42). In 
the region of interest say one chord upstream and two chord 
downstream of the stator, let us assume the amplitude U 

wr .O) 

in Eq. (42) to be constant and equal to its value at the stator 
midchord station. Such an approximation would yield higher 
estimates but is convenient for comparison with the correspon- 
ding "ordered" dipole radiation from periodic fluctuating forces 
on the outlet guide vanes discussed in section Il.C.l. 


Considering the radiation from the "ordered" quadrupoles in the 
region upstream of the stator, we can write the rms value of 
acoustic pressure as 


OO 


4P 0 ) = 575 p “wr K|) cos 1 '*' J B -v< krsin 'J' > V 


P s Hi) 


dx' 


(55) 


Similarly from the region downstream of the stator we can write 

2c_ 


“ 572 J B . v (krsin^)rdr 


f 8 

J ws i(i> 


dx" 


(56) 


o 

From Eq. (23) describing the acoustic radiation from "ordered" 
dipoles on the outlet guide vanes due to rotor wake impingement. 
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we can restrict our attention to only the fundamental of the 
axial component of the fluctuating load, and write 

— IfTT 

^ P (t) = 4iR72 (4L) (0 SlnX s COS ^ J B _ v (krsin ^') < 57 > 

/\ 


where ( 4L)^ is the amplitude of the fundamental of fluctu- 
ating lift force on the vane element in the annulus 2urdr 
due to periodic impingement by the rotor wakes. 

The magnitude of depends upon the amplitude of the 

fundamental of the upwash (component normal to \p fluctu- 
ations experienced by the vane element, which can be approxi- 
mated by 

sin( v+v) A 

i — 2_ u 

cos wr,„, 


and the lift transfer function discussed by Kemp and sears in 
ref. (11). 


Using the value suggested in ref. (19) for the lift transfer 
function we can write 


/n a sin( v _+ v) 

( 4l) = U c dr — — 

(i) P o wr s cos v 

' 1 0) s 


1 + 


Trka 0 c 

V 

o 


rh 


(58) 


Let us compare the contribution from the ordered quadrupoles 
considered in Eq. (55) with that from the ordered dipoles as 
given by Eq. (57). Substituting Eq. (58) into (57) we can 
obtain the ratio 
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AP {]) from Eq. (55) 


COS V 

s 


„ k r , f , 

= 4 — — cos^ — . , . , „ . 

ip,, from Eq. (57) v c s sln * B Sln( VV 


( 1 ) 


X . 


1 + 


'TTka 0 c g 


V o 


00 

J< /• A 

/ J 


10 ) 


V 

o 


dx* 


( 59 ) 


Let us assume U 


ps 




at the leading edge of the stator to be 


approximately equal to 0.4^ and obtain 


OO 

/ 


A 

U 


P S KD . , n . a s 
V dx = 0,4 2 tt 


We can also use simplifying approximations such as 
( + v) ^ it/2 

x S 


cos « 


sin\ tan» 


U 


U. 


For conditions usually met in lifting fans we can write 
f ka 0°s| 

I 1 *"— } 

Substituting the above approximations into Eq. (59) we obtain 
the ratio 


4P (I ) from Eq. ( 55 ) 
from Eq.O?) 


U 

, B r 0 ’ 

= 1.6 y M t cr s cos Mf — 


(60) 


and observe that the contribution from the quadrupoles considered 
here is of the same’ order of magnitude as that from the periodic 
forces on the vanes due to impingement by rotor wakes. 
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By comparing the contributions from the ordered quadrupoles 
in the downstream and upstream region of the stator as given 
by Eqs. (56) and (55) respectively, we obtain a ratio similar 
to that derived in earlier Eq. (52). Consequently, we can 
conclude that the ordered quadrupoles in these two regions are 
of equal importance in evaluating fan noise. 

III.C. Random Quadrupoles Due to Inflow Turbulence 
In the preceding subsection we have discussed the acoustic 
radiation from only the velocity fluctuations which are periodic. 
Since velocity fluctuations can also contain a random component 
due to turbulence, one should consider the resulting random 
quadrupole noise sources. For the present let us restrict our 
attention to the turbulent fluctuation in the inflow to the fan. 
The random fluctuations by themselves give rise to a broad 
band noise, whose magnitude is governed by the intensity and 
power spectrum of the inflow turbulence. However, the presence 
of the periodic component related to the rotor flow field can 
give rise to noise at blade passing frequency and its harmonics. 

III.C.l. Effect of Rotor Potential Flow Field 
In the upstream region of the rotor, we have periodic velocity 
fluctuations arising from potential flow around the blades, as 
described in subsection III. A. 2, and random velocity velocity 

fluctuations whose power spectrum is discussed in subsection 
III. A. 4. Let us consider a point (r,0,£) in- the fan annulus 
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and an elemental volume rdrd0d£ within which the velocity can 
he considered to be the same as at the point. Let us restrict 
our attention to the axial components of velocity perturbation 
and substitute 

° = U pr, + v °, 

into Eq. (25) to define the elemental quadrupole strength. 
Considering the power spectrum 


<t>(±oj) = v? — 
' ' 0 ITU 




for the turbulent fluctuations, and the sinusoidal nature 
U = G cosB£>t 

P 1(1) P i(i) 

of the rotor related component at the blade passing frequency 
we can show that the power spectrum of the elemental quadrupole 
source is 


A 


V ±U) = 2pU pr 


v 


■Ml 0 




(61) 


The convection velocity U c of the "frozen convective pattern" 
model of turbulence employed in the above equation can be 
approximated by the mean axial velocity U . Let us consider 

°i 

an eddy size of dimensions £ , X 2 and along the axial. 


- 43 - 



radial and circumferential directions respectively. For iso- 
tropic turbulence we have £ = 2X 7 = 2£ 3 = £. A simplifying 
assumption would be that the turbulent fluctuations are perfectly 
correlated within the eddy volume and uncorrelated with fluctu- 
ations elsewhere. However, the rotor related sinusoidal fluc- 
tuations give rise to different retarded times for radiation 
from different portions of the eddy. Acoustic radiation from 
such quadrupole sources was discussed in ref. (20). The mean 
square pressure at the field point due to sources with an 
annulus 2TTrdr extending from rotor leading edge to infinite 
distance upstream can be written as 


cos & v 


ITUn 


AP 2 = (ka o r £ ] £ 2 £ 3 4> q (B£) f,f 3 2ir r dr 


327T a^ 


A 


hA 


where f, = / {e 


-mr( 2x'+A )/d 


-i 


dx' 


+ J e 

4a 

J d . tt^ 3 
= ~» — 


-4rrx*/d r rt£ f /d -TT^/d 


f 


-e 




dx' 


and <f> (B Q) = the right hand side of Eq, (61) at co= B Q 


By integrating above equation from r = r h to r fc one can 
evaluate the acoustic radiation from longitudinal-axial 
quadrupole sources in the fan annulus upstream of the rotor. 
Results of such computations on Rao Fan 2 at 5 ft. distance 
on axis are presented in Fig. (17) for various values of 
length scale £. 
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III.C.2. Effect of Rotor Wakes 

Downstream of the rotor blades, we have periodic fluctuations 
in velocity as discussed in section III. A. 3. Similar to the 
analysis presented in the preceding subsection, let us examine 
the effect of the axial component of the fundamental of such 
velocity fluctuations. Up to an axial downstream distance of 
twice the rotor chord, the axial component of the fundamental 
of velocity fluctuations have already been described by Eq. (42). 
The periodic component of the velocity fluctuations considered 
in the preceding section are those described by Eq. (41). We 
recognize the similarity between these two equations and the 
acoustic radiation from the random quadrupole source in this 
region can be approximated by comparing the integral 

oo 



with the integral 



(63) 


(64) 


where is the velocity relative to the rotor blades. 


We note that the ratio of the above two integrals is approx- 
imately unity when we consider the values of C and a 

D r u r 

usually occurring in axial fans. Consequently, we can estimate 
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that the noise radiated from the random quadrupole sources 
present in this downstream region is of the same order of 
magnitude as that discussed in the preceding subsection. 

III.C.3. Effect of Stator Flow Field 

The perturbations in flow field caused by the presence of stator 
vanes are of a steady-state type even though they depend upon 
spatial location. Consequently, the effect of random fluctu- 
ations associated with turbulence will give rise to quadrupole 
radiation having the same spectral distribution as the turbulence 
component. Our attention presently being the n6ise at blade 
passage frequency and its harmonics, consideration of the effects 
referred above are beyond the scope of this report. 

III.D. Random Quadrupoles Due to Turbulence in the Wake Regions 
In addition to the homogeneous turbulence present in the inflow 
to the fan, we have random turbulent fluctuation of velocity in 
the wake regions of the various elements of the fan. Even though 
the intensity of turbulence in these regions is greater than 
that in the inflow it is present only in a small portion of the 
fan annulus. Because of the rotor related periodic fluctuations, 
the associated quadrupoles can generate noise at blade passage 
frequency and its higher harmonics. Evaluation of acoustic 
radiation from the wake regions will be included in our future 
investigations. 
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IV. CONCLUDING REMARKS 

This interim report describes the analytical methods employed 
along with some preliminary results obtained in our effort of 
evaluating noise sources in axial fans. Since acoustic meas- 
urements on fans show prominent peaks at blade passing fre- 
quency and its harmonics ; we directed our attention to the 
discrete frequency noise. The information presented in this 
report is a' step towards identifying predominant noise sources 
in a fan, and must be regarded as far from a complete under- 
standing of the fan noise problem. However, based on the 
representative results given in this report the following con- 
clusions can be reached. 

(1) The effect of steady-state blade loading can be 
ignored in estimating noise from fans with multi- 
bladed rotors. 

(2) The periodically fluctuating forces on rotor blades 
or stator vanes arising from wake-impingement can be 
a dominant noise source depending upon the blade- 
vane ratio. The noise level as well as its direc- 
tivity can be altered by changing the number of 
stator vanes. 

(3) Quadrupole radiation from periodic fluctuations of 
axial velocity components in rotor-stator inter- 
action can be of the same order of magnitude as the 
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dipole radiation from the periodic fluctuating loads 
on blades or vanes caused by wake impingement. 

(4) Inherent cancellation of acoustic radiation from 
"ordered" sources distributed in the fan annulus 
yields zero noise along fan axis. 

(5) Turbulence present in the inflow to the fan can be 
an important contributor to discrete frequency noise 
because of the random forces on the rotor blades. 
Acoustic radiation from random quadrupoles is rela- 
tively unimportant. 

As mentioned earlier, further effort on several items briefly 
mentioned in this report is necessary for a complete under- 
standing of acoustic radiation from axial fans and to be able 
to arrive at design features to reduce noise. 
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Fig. 1. Field Point Location Relative to the Source Point 
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Fig. 4 


Sound Pressure Levels at Bla< 
State Blade Loading in R 
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Sound Pressure Levels at Blade Passing Frequency Computed from Steady 
Blade Loading for LF336 (150 ft. distance from fan). 



Fig. 5 





phase 



Fig. 7. Radial Variation of the Phase Angle at Fundamental (a=l) of 
the Fluctuating Lift Force on Blade Element in Rao Fan 2. 
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Sound Pressure Levels at Blade Passing Frequency Computed from 
Periodic Blade Loading in Rao Fan 2 (5 ft, ditance from fan) , 


Fig. 8 
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Fig. 9. Radial Variation of ath Component Amplitude 
of Fluctuating Lift Force on Blade Element in 
LF336 Due to Major Strut 



phase angle 



Fig. 10. Radial Variation of a- 
of Fluctuating Lift F< 
LF336 Due to 






Fig. 11. Vectorial Representation of the ath Component 
of Fluctuating Lift Force on Blade Element in 
LF336 Due to Major Strut 


aspect angle xjf 


Fig. 12. Sound Pressure Levels at Blade Passing Frequency Computed from Periodic 

Blade Loading in LF336 at Design Speed 
(effect of major strut wake only) 
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Fig. 14. Acoustic Radiation on 


( inlet turbulence 


frequency kHz 

the Axis from Random Forces on Rotor Blades 
in LF336 

tensity = 3 %. distance from fan = 150 ft.) 







(ZlL) lbs. per ft. length 








SPLdB ref. 0.0002 U bar 



Fig . 18 . Sound Pressure Levels at Blade Passing Frequency Computed from 
Periodic Vane Loading in LF336 (150 ft. distance from fan). 
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